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Background: Heavy-flavor production in p+p collisions is a good test of perturbative-quantum- 
chromodynamics (pQCD) calculations. Modification of heavy-flavor production in heavy-ion colli- 
sions relative to binary-collision scaling from p+p results, quantified with the nuclear-modification 
factor (Raa), provides information on both cold- and hot-nuclear-matter effects. Midrapidity heavy- 
flavor Raa measurements at RHIC have challenged parton-energy-loss models and resulted in upper 
limits on the viscosity /entropy ratio that are near the quantum lower bound. Such measurements 
have not been made in the forward-rapidity region. 

Purpose: Determine transverse-momentum, pr spectra and the corresponding Raa for muons 
from heavy-flavor mesons decay in p+p and Cu+Cu collisions at ^/snn = 200 GeV and y — 1.65. 

Method: Results are obtained using the semi-leptonic decay of heavy-flavor mesons into negative 
muons. The PHENIX muon-arm spectrometers measure the pr spectra of inclusive muon candidates. 
Backgrounds, primarily due to light hadrons, are determined with a Monte-Carlo calculation using a 
set of input hadron distributions tuned to match measured-hadron distributions in the same detector 
and statistically subtracted. 

Results: The charm-production cross section in p+p collisions at i/s = 200 GeV, inte- 
grated over pr and in the rapidity range 1.4 < y < 1.9 is found to be da C c/dy — 0.139 ± 
0.029 (stat) 1q Q58 ( s y s t) m b- This result is consistent with a perturbative fixed-order-plus-next-to- 
leading-log (FONLL) calculation within scale uncertainties and is also consistent with expectations 
based on the corresponding midrapidity charm-production cross section measured by PHENIX. 
The Raa for heavy-flavor muons in Cu+Cu collisions is measured in three centrality intervals for 
1 < pr < 4 GeV/c. Suppression relative to binary-collision scaling (Raa < 1) increases with 
centrality. 

Conclusions: Within experimental and theoretical uncertainties, the measured heavy-flavor 
yield in p+p collisions is consistent with state-of-the-art pQCD calculations. Suppression in cen- 
tral Cu+Cu collisions suggests the presence of significant cold-nuclear-matter effects and final-state 
energy loss. 

PACS numbers: 25.75. Dw 



I. INTRODUCTION 

Understanding the energy loss mechanism for par- 
tons moving through the hot dense partonic matter pro- 
duced in heavy-ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) and the Large Hadron Collider 
(LHC) is a key priority in the field of heavy-ion colli- 
sion physics p], . Production of heavy quarks in heavy- 
ion collisions can serve as an important tool for better 
understanding properties of the dense matter created 
in such collisions. In particular, because of their large 
mass, heavy quarks are almost exclusively produced in 
the early stages of heavy-ion collisions and can therefore 
serve as a probe of the subsequently created medium. 
The large mass scale associated with the production of 
heavy quarks also allows one to test perturbative Quan- 
tum Chromodynamics (pQCD) based theoretical models 
describing high energy collisions. 

Recent measurements of heavy-quark production in 
heavy-ion collisions [H-Q exhibit a suppression, which 
is larger than expected and not easily reconciled with 
early theoretical predictions [U, 0] • In these calculations 
the dominant energy loss mechanism for fast partons is 
gluon bremsstrahlung [8, 9]. In this context, it was pre- 
dicted that heavy quarks would lose less energy than 
light quarks due to the so-called dead-cone effect 
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The disagreement between this prediction and experi- 
mental results led to a consideration of alternative in- 
medium parton energy loss mechanisms, assumed earlier 
to have a small effect on heavy quarks compared to ra- 
diative energy loss. In particular, it was suggested that 
heavy quarks can lose a significant amount of their en- 
ergy through elastic collisions with in-medium partons 
(collisional energy loss mechanism) (Tll - H^ . especially 
in the intermediate transverse momentum range (pt ~ 

3 8GeV/c) in which most of the RHIC open heavy 

flavor measurements are performed. Additional mecha- 
nisms for in-medium e nerg y loss for heavy quarks have 
also been suggested [i~4lll5||. Despite recent progress, still 
needed is a universal theoretical framework describing 
precisely the production of heavy quarks and their sub- 
sequent interactions with the partonic medium created in 
heavy-ion collisions. Also needed are accurate measure- 
ments of heavy-quark production in heavy-ion collisions, 
which are critical to test and constrain the theoretical 
predictions. 

Hidden- heavy- flavor (J/ip) production has also been 
extensively measured in heavy- ion collisions [l6l.[r7j. The 
production of J/tp mesons is expected to be affected by 
the formation of a quark-gluon plasma due to the in- 
terplay of several competing mechanisms, including sup- 
pression due to a color scree ning mechanism similar to 
the Debye screening in QED [18] and enhancement due 
to the coalescence of uncorrelated cc pairs from the hot 
medium [19h21| . The magnitude of such an enhance- 
ment depends strongly on the production cross section of 
open-heavy flavor in heavy-ion collisions, measurements 
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of which are therefore essential to the interpretation of 
heavy quarkonia results. 

A well-established observable for quantifying medium 
effects in heavy-ion collisions is the nuclear-modification 
factor, Raa- 



Raa 



1 Q-aa 

-Wcoll CTpp 



(1) 



where ctaa and the invariant cross sections for 

a given process in A + A collisions and p + p collisions, 
respectively, and iV co u is the average number of nucleon- 
nucleon collisions in the A + A collision, evaluated using 
a simple geometrical description of the A nucleus (22[. 

For processes that are sufficiently hard (characterized 
by large energy transfer), Raa is expected to be equal to 
unity in the absence of nuclear effects. A value smaller 
(larger) than unity indicates suppression (enhancement) 
of the observed yield in A + A collisions relative to ex- 
pectations based on p + p collision results and binary- 
collision scaling. 

Open-heavy-flavor production has been measured by 
the PHENIX experiment at midrapidity (\n\ < 0.35) 0. 
This paper presents the measurement of open-heavy- 
flavor production at forward rapidity (1.4 < \rj\ < 1.9) 
in Cu + Cu and p + p collisions, and the resulting Raa, 
using negatively-charged muons from the semi-leptonic 
decay of open-heavy- flavor mesons. 

The paper is organized as follows: Section|TT]prcscnts a 
short overview of the PHENIX detector subsystems rel- 
evant to these measurements followed by a description 
of the data sets and track selection criteria. Section IIIII 
presents a detailed description of the methodology for 
measuring the invariant cross section in p + p collisions 
and Raa in Cu + Cu collisions for muons from heavy- 
flavor- meson decays. Results are presented in Section Hvl 
and compared to existing measurements as well as theo- 
retical predictions in Section IVl 



ZDC South 




Side View 



FIG. 1: (color online) Side view of the PHENIX muon detec- 
tors (2005). 



in the MuTr to the penetration depth of the particle in 
the MuID (that is, the last MuID gap a given particle 
reaches) is the primary tool used to identify muons with 
respect to the residual hadronic background. Measured 
muons must penetrate 8 to 11 interaction lengths in total 
to reach the last gap of the MuID. This corresponds to 
a reduction of the muon longitudinal momentum (along 
the beam axis) of Sp z = 2.3(2.45) GeV/c in the south 
(north) muon arm. The MuID is also used in the online 
data acquisition to trigger on collisions that contain one 
or more muon candidates. 

Also used in this analysis are the Beam-Beam Coun- 
ters (BBC) [24J, which comprise two arrays of 64 quartz 
Cerenkov detectors that surround the beam, one on each 
side of the interaction point. The BBCs measure charged 
particles produced during the collision in the pseudora- 
pidity range 3 < \rj\ < 3.9 and determine the collision's 
start-time, vertex longitudinal position, and centrality 
(in Cu + Cu collisions). The BBCs also provide the min- 
imum bias trigger. 



II. EXPERIMENTAL SETUP AND DATA SETS 

A. The PHENIX Experiment 

The PHENIX experiment is equipped with two muon 
spectrometers [23[, shown in Fig. [TJ that provide pion 
rejection at the level of 2.5 x 10~ 4 in the pseudora- 
pidity range —1.2 < rj < —2.2 (south muon arm) and 
1.2 < r) < 2.4 (north muon arm) over the full azimuth. 
Each muon arm is located behind a thick copper and iron 
absorber and comprises three stations of cathode strip 
chambers (the Muon Tracker, or MuTr), surrounded by 
a radial magnetic field, and five "gaps" (numbered 0- 
4) consisting of a plane of steel absorber and a plane 
of Iarrocci tubes (the Muon Identifier, or MuID). The 
MuTr measures the momentum of charged particles by 
tracking their motion in the surrounding magnetic field. 
Matching the momentum of the particles reconstructed 



B. Data Sets 

Two data sets, recorded in 2005, are used in this 
analysis: p + p collisions and Cu + Cu collisions at a 
center of mass energy per nucleon-nucleon collision of 
= 200 GeV. 

The p + p data used for this analysis have been 
recorded using two muon enriched triggers, in coinci- 
dence with the Minimum Bias (MB) trigger, which re- 
quires at least one hit in each of the BBCs and covers 
approximately 55% of the total p + p inelastic cross sec- 
tion. These two muon triggers rely on the information 
recorded in the MuID. The first (Deep) trigger requires 
one or more muon candidates to reach the last plane of 
the MuID (Gap4), whereas the second, less strict, (Shal- 
low) trigger requires one or more muon candidates to 
reach at least the third MuID gap (Gap2). The inte- 
grated luminosity sampled with these triggers and used 
for this analysis is 44.3 nb _1 (48.7 nb _1 ) for the south 
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TABLE I: Centrality characterization variables for Cu + Cu 
collisions. 



centrality 


0-20% 


20-40% 


40-94% 




151.8 ± 17.1 


61.6 ± 6.6 


11.23 ± 1.3 


-/Vpart 


85.9 ± 2.3 


45.2 ± 1.7 


11.7 ± 0.6 



(north) muon arm. 

All Cu + Cu data used for this analysis have been 
recorded using the Minimum Bias trigger described 
above. For Cu + Cu collisions, this trigger covers ap- 
proximately 94 % of the total inelastic cross section. The 
integrated luminosity sampled with this trigger and used 
for this analysis is 0.13 nb _1 , using a total Cu + Cu in- 
elastic cross section seen by the minimum bias trigger 



el 

'Cu+Cu 



= 2.91 b. 



C. Centrality Determination 

The centrality of each Cu + Cu collision is determined 
by the number of hits in the BBCs. Three centrality 

bins are used for this analysis: 20%, 20 40% 

and 40 94%, where 20% represents the most 

central 20% of the collisions. For a given centrality, the 
average number of nucleon-nucleon collisions (N co \\) and 
the average number of participating nucleons (Af part ) are 
estimated using a Glauber calculation [22j coupled to a 
model of the BBC response. Values of N co \\ and N pal - t 
for the three centrality bins defined above are listed in 
Table H 

To ensure that the centrality categories are well de- 
fined, collisions are required to be within ±30 cm of the 
center of the PHENIX detector along the beam axis. 



Z=0 plane 



Ay 



(x,y) = (0,0) 




Track associated MuID road 
projection to Z=0 plane 



FIG. 2: (color online) Schematic representation of r re f vari- 
able. 
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FIG. 3: (color online) Schematic representation of track se- 
lection variables DG0 and DDG0. 



D. Track Selection 

This section outlines the track-based selection vari- 
ables. 

zbbc The event vertex longitudinal position is measured 
by the BBC detector. For low-momentum tracks 
(pt < 2GeV/c) reconstructed in north (south) 
muon arm we demand zbbc > (zbbc < 0). This 
arm-dependent cut improves the signal to back- 
ground ratio because light hadrons produced dur- 
ing the collision have a probability to decay into 
a muon that increases with their distance from 
the front muon arm absorber, whereas muons from 
short-lived heavy-flavor hadrons have a yield that 
is independent of zbbc ( see al so Section Ull C 3|) . 

Zflt The vertex longitudinal position of a track evaluated 
using a fit of the track position and momentum 
measured in the MuTr and extrapolated backward 



through the front absorber towards the interaction 
point, together with the BBC vertex measurement. 

Ar hfts Tr Tnc total number of track hits in the three MuTr 
stations. A given track can have up to 16 MuTr 
hits. 

N hiti ID Tnc total number of track hits in the five MuID 
gaps. A given track can have up to 2 hits in each 
gap (10 in total). 

r re f The distance to the beam axis of the track, as re- 
constructed in the MuID only, when extrapolated 
(backward) to z = (illustrated in Fig. [2]). 

Road Slope The slope of the track, as reconstructed 
in the MuID only, measured at MuID GapO: 
\J (dx/dz) 2 + (dy/dz) 2 . A cut applied to this vari- 
able eliminates combinatorial background gener- 
ated in the high hit-density region closest to the 
beam pipe. 
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DGO The distance between the track positions calcu- 
lated in the MuTr and in the MuID, evaluated at 
the MuID GapO z position (see Fig. 03). 

DDGO The difference between the track angles calcu- 
lated in the MuTr and in the MuID, evaluated at 
the MuID GapO z position (see Fig. 02). 

5z The difference between the event vertex longitudinal 
position reconstructed by the BBC (zbbc) and the 
track longitudinal position provided by the track 
reconstruction algorithm: Sz = |zbbc — Zfitl- 

p59 the effective scattering angle of the track in the front 
absorber, 59, scaled by the average of the momen- 
tum measured at the vertex and at MuTr Station 
1: p = (pvtx +Psti)/2, where 59 is given by: 

SO = cos-i ( l^lstA (2) 

V Pvtx-Pstl / 

where Psti ^ s the momentum vector measured at 
Station 1 and p vtx is the momentum vector at 
the vertex. For a given track, 56 essentially mea- 
sures the track deflection in the front absorber due 
mostly to multiple scattering and radiative energy 
loss, but also to the magnetic field upstream of sta- 
tion 1. This deflection is expected to be inversely 
proportional to the track total momentum. Scal- 
ing the scattering angle 69 by the track momentum 
therefore ensures that the p59 distribution is ap- 
proximately Gaussian with a constant width for all 
Pt bins. 

Cut values applied to these variables are, in some 
cases, pt-, species- and/or centrality-dependent. Within 
a given pt, species and centrality bin, the same cut val- 
ues are applied to both Monte Carlo simulations and real 
data. 

Even after all cuts are applied to select good quality 
muon candidates, there remains a small contamination 
of misreconstructed tracks caused by: 

• Accidental combinations of hits in the muon tracker 
that do not correspond to a real particle. 

• Tracks arising from interactions between the beam 
and residual gas in the beam pipe or between the 
beam and beamline components. 

These misreconstructed tracks, later denoted Np, are 
not completely reproduced by experimental simulations 
and must be estimated and properly subtracted from the 
inclusive muon sample to evaluate the amount of muons 
from heavy-flavor decay. The method by which Np is 
estimated is based on the distributions of the p59 and 5z 
variables and is described in more detail in Section [ill Bl 

Note: positive muons are not used in this analysis due 
to a poorer signal/background ratio resulting from the 
fact that both anti-protons and negative kaons are more 



strongly suppressed by the MuTr front absorbers than 
their positive counterparts. The rapidity interval used for 
this measurement is smaller than the rapidity coverage 
of the PHENIX muon spectrometers (1.2 < \r}\ < 2.2) to 
reduce uncertainties in the acceptance calculation. 

III. METHOD FOR THE MEASUREMENT OF 
HEAVY-FLAVOR MUONS 

A. Overview 

The methodology used to measure heavy-flavor muon 
(i.e., muons from heavy-flavor meson decay) production 
in p + p and Cu + Cu collisions is described in this sec- 
tion. This analysis is a refinement of techniques originally 
developed in (25^271 . 

For both p + p and Cu + Cu collisions the double dif- 
ferential heavy flavor muon invariant yield is defined by: 

d 2 N» = 1 Ni - N c - N F 

2Trp T dp T dr] 27rp T Ap T A?7 N evt e%xc Ae 

where Nj is the total number of muon candidates in 
the bin, consisting of the tracks that reach the last gap 
of the MuID (Gap4) and pass all track selection crite- 
ria; Np is the estimated number of misreconstructed 
tracks that pass the track selection cuts accidentally (Sec- 
tion IIIIBI) ; Nq is the number of tracks corresponding to 
the irreducible hadronic background, as determined us- 
ing a hadron cocktail approach (Section l"!!! C|) : iV cvt is the 
number of events, Ae is the detector acceptance and ef- 
ficiency correction (Section [III El) , and Cg^f is the BBC 
trigger efficiency for events in which a heavy-flavor muon 
at forward rapidity is present. This efficiency amounts 
to 79% (100%) in p + p (Cu + Cu) collisions. 

The p + p and Cu + Cu invariant yields determined 
with Eq. 03 can be used directly to determine the heavy- 
flavor muon i?AA (Eq. EQ)- However, in order to mini- 
mize the systematic uncertainty associated with the es- 
timate of the hadronic background by canceling the part 
of this uncertainty that is correlated between the p + p 
and the Cu + Cu analyses, Raa is calculated separately 
for a given i th version of the Monte-Carlo simulation of 
hadron cocktail used in the estimate of Nc'- 

R >- aa= 1 f d 2 N C n+Cu/dp T dri Y 
AA Ncoii V d 2 N p+p /dp T dr] J 

The final value for Raa is then determined by taking the 
mean of the values obtained for the different cocktails, 
each weighted by its ability to reproduce measured data, 
as discussed in Section UlI CI 

B. Contamination from Misreconstructed Tracks 

Np, the number of misreconstructed tracks that ac- 
cidentally pass all track quality cuts, is estimated using 
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3.0 GeV/c < p T < 4.0 GeV/c 

■ Raw pS9 distribution, |5z|<2cm 
▼ Raw pSe distribution, |Sz|>2cm 



3.0 GeV/c < p T < 4.0 GeV/c 

■ Raw p86 distribution, |5z|<2cm 
▼ Scaled p59 distribution, |5z|>2cm 



FIG. 4: (color online) pSQ distributions for north arm inclu- 
sive muon candidates, 3 < pp < 4GeV/c. The top panel 
compares the distribution inside (black squares) and outside 
(red triangles) the Sz cut. The bottom panel compares the 
same distributions, but the distribution outside the Sz cut 
(red triangles) is normalized to the distribution inside the Sz 
cut (black squares) in the region pS8 > p5(9 max - In both pan- 
els, the vertical dashed line corresponds to p5f? max . 



Using Equations [5] and |6l it is found that Np amounts 
to less than 1% of the inclusive muon sample in the lowest 
pp bin (1 < pp < 1.5 GeV/c) and increases with pp up 
to about 5% for the highest pp bins. Uncertainties on 
these estimations are negligible in the final results. 



C. Hadron Cocktail 

Charged pions and kaons are the largest source of par- 
ticles in the PHENIX muon arms. Other species (p, 
p, K®, Kl) have small but nonzero contributions. Al- 
together, these light hadrons constitute the main back- 
ground source for the measurement of muons from heavy- 
flavor meson decay. 

One can define three contributions to this background, 
depending on how the particles enter the muon spectrom- 
eter: 



Decay muons - light hadrons that decay into muons 
before reaching the first absorber material. Since 
these particles enter the spectrometer as muons, 
a fraction of them also penetrate all the absorber 
layers of the MuID and enter the pool of inclusive 
muon candidates. 

Punch-through hadrons - hadrons produced at the 
collision vertex that do not decay, but penetrate all 
MuID absorber layers, thus also being (incorrectly) 
identified as muons. 



the p56 distribution inside and outside of the Sz cut de- 
fined in Section III Dl These two distributions are shown 
in the top panel of Fig. 2) The distribution inside the 
Sz cut (black squares) shows two contributions: a peak 
at pS8 = 0.05rad-GeV/c, corresponding to the expected 
multiple scattering of muons in the front absorber, and 
a tail out to large values of pSd. In the distribution out- 
side the Sz cut (red triangles), the signal peak has dis- 
appeared, and only the tail remains. Note that the tail 
extends below the pSQ cut; this is the Np contribution. 
Using the fact that the shape of this tail appears to be 
the same on both sides of the Sz cut, one can estimate 
Np using: 

N F = aN F ' (5) 

where Np' is the number of tracks with pS9 < p59 max 
but Sz > Sz max , and a normalizes the tails of the two 
distributions above the pS9 cut: 

= N(pS9 > p59 majK ,Sz < fe max ) 
N(j>56>p68 max ,6z>6z max ! 

The bottom panel of Fig. 0] shows the p88 distribution 
inside the Sz cut (black squares, identical to the corre- 
sponding distribution in the top panel) and the distribu- 
tion outside the Sz cut (red triangles from the top panel) 
after scaling by a (Eq. [B]). 



Decay-in-MuTr - hadrons produced at the collision 
vertex that penetrate the muon arm front absorber 
and decay into a muon inside the MuTr tracking 
volume, with the decay muon then passing through 
the rest of the MuTr and the MuID. Most such par- 
ticles are simply not reconstructed because of the 
decay angle between the primary hadron and the 
decay muon. However, some can be reconstructed, 
usually with an incorrect momentum assigned to 
the track. Due to the exponential pp distribution, 
even a small number of such tracks can form a 
significant background at high pp, but for the pp 
range in this analysis this contribution is small. 

While decay muons can not be distinguished from 
punch-through hadrons and heavy-flavor muons on an 
event- by- event basis, their production exhibits a strong 
vertex dependence, as illustrated in Fig. [5] This feature 
plays a key role in constraining heavy-flavor background 
(Section MHM- 

A series of Monte Carlo simulations ( "hadron cocktail 
packages") are used to estimate the overall background 
due to light hadron sources. The construction of a given 
hadronic cocktail package involves the following steps: 

1. Generate a primary hadron sample based on 
parameterized pp and y distributions (Sec- 
tion men)- 
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FIG. 5: (color online) Vertex z distribution of muon candi- 
dates reconstructed in north (z > 0) MuID Gap4, relative 
to the event vertex z distribution (black circles). The vertex 
z dependencies of the various contributions to the inclusive 
muon spectra are represented schematically as colored boxes. 



2. Propagate these hadrons through the muon spec- 
trometer using the complete GEANT3 28] PHENIX 
simulation. Each hadron cocktail package uses 
one of the two hadron shower codes provided by 
GEANT3: G-fluka or GHEISHA with a scaled value 
of the hadron-Iron interaction cross section (Sec- 
tion nrcii). 

3. For the Cu + Cu analysis the simulated hadrons are 
then embedded in real events in order to account 
for deterioration of the reconstructed track quality 



due to high hit multiplicity. 

4. Reconstruct the resulting particles using the same 
reconstruction code and track quality cuts used in 
the real data analysis. ( Section III D|) . 

5. Tune (that is, re-weight) the input px distributions 
(from step 1) to match hadron distributions mea- 
sured in the muon arm (Section 1111 C 31) . 



1. Input Particle Distributions 

Particle distributions required as input to the hadron 
cocktail have not been measured over the required y and 
Pt range at RHIC energies. We therefore use a combi- 
nation of data from PHENIX, BRAHMS and STAR, to- 
gether with Next-to-Leading Order (NLO) pQCD calcu- 
lations to derive realistic parameterizations of these dis- 
tributions. An exact match to actual distributions is not 
necessary since the input distributions are re- weighted to 
match measured hadron distributions before being used 
to generate estimates of Nc (Section [ill C 3|) . 

We start with the 7r° spectrum in p + p collisions at 
y = measured by PHENIX [H| . This is extrapolated to 
y = 1.65 in two steps. First, an overall scale factor is ob- 
tained from a Gaussian parameterization of the charged 
pion dN/dy distribution measured by BRAHMS [30|. 
Next, the px shape is softened using a parameterization 
of the ratio of unidentified hadron pt spectra measured 
by BRAHMS at rj = and 77 = 1.65 [HI, HI]: 



dN/dp T (TT ± ,y= lM) = dN/dp T {n°,y = 0) x exp(-i(1.65/2.25) 2 ) x (1 - (0.1 ■ p T [GeV/c] - 1)) (7) 

I 



Next we extrapolate this spectrum over the range 
1-0 < y < 2.4 using a series of Next-to-Leading 
Order (NLO) calculations [33j to obtain the ratio 
dN/dpT^ ± ,y)/dN/dp T (^ ± ,y = 1-65). Figure H shows 
a comparison of the hadron cocktail input for charged 
pions compared to charged-pion distributions at y = 
and y = 2.95. Spectra for other hadron species in the 
cocktail are obtained by multiplying the parameterized 
pion spectra by parameterizations of measured values of 
hadron-to-pion ratios, as a function of px- 

With 8-11 interaction length of material prior to MuID 
Gap4, approximately 4000 hadrons must be simulated to 
obtain a single hadron reconstructed as a muon. Given 
this level of rejection, it is very CPU intensive to generate 
a sufficient sample of high px hadrons using realistic pr 
spectra. A standard technique is to throw particles with 
a flat pt spectrum and then weight them with a realis- 
tic distribution. However, interactions in the absorber 



in front of the MuTr and decays in the MuTr volume 
can both result in particles being reconstructed with in- 
correct momentum. Due to the steeply falling nature of 
the pt spectrum, tracks with low momentum and incor- 
rectly reconstructed with a higher momentum can have a 
significant contribution at high pt, with respect to prop- 
erly reconstructed tracks. As a compromise designed to 
ensure statistically robust samples of both tracks with 
initial high pt and with misreconstructed high pt, we 
multiply the realistic pt distributions by pt 2 to form 
the simulation input pt distributions, and re-weight the 
output of the simulation by 1 /pt 2 to recover the initial 
distribution. 

The particles in the primary hadron sample used as 
input to each hadron cocktail package are generated as 
follows: 

• The particle type and rapidity are chosen based 
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FIG. 6: (color online) Pion cross sections as a function of 
Pt used as initial hadron cocktail input, for several rapid- 
ity intervals in [1.0, 2.2] (blue lines) compared to a fit to the 
PHENIX 7T° data at y = [H (black line, open black circles) 
and BRAHMS tt~ data at y = 2.95 [HI (open black circles). 

on dN/dy values obtained by integrating the un- 
weighted pt distributions described above. 

• The particle's transverse momentum is chosen 
within the range 0.8 < pt < 8 GeV/c using the 
p T 2 _weighted pt distributions described above. 

• Since the muon spectrometer acceptance shows lit- 
tle dependence on the vertex z position, the parti- 
cle's z origin is chosen from a flat distribution over 
the range — 35 < z < 35 cm. 

• The particle's azimuthal angle, <p, is chosen from a 
flat distribution over 2-k. 



2. Hadron Cocktail Packages 

Modeling hadron propagation through thick material 
is known to be difficult and neither hadron shower code 
available in geant3 (g-fluka and gheisha) is able to 
reproduce measured data in the PHENIX muon arms. 
The approach we have chosen to circumvent this issue 
is to produce a range of background estimates using a 
set of hadron cocktails (referred to as packages) , each of 
which uses one of the GEANT hadron shower codes and 
a different, modified, value of the hadron-Iron interac- 
tion cross section. The set of background estimates are 
then combined in a weighted fashion to extract central 
values for production yields, Raa, and the contribution 
to the systematic uncertainty on these quantities due to 
the uncertainty in hadron propagation. 

Using the default hadron-ion cross section, fluka sim- 
ulations produce more muon candidates than gheisha 



simulations, therefore the FLUKA cross sections are in- 
creased relative to the default and the gheisha cross sec- 
tions are decreased. The cross section modifications are 
referred to in terms of percentage, so that a 6% increase 
is referred to as 106%. Five packages are used in this 

analysis: FLUKA105 (or FL105), FL106, FL107, GHEISHA91 

(or GH9l) and GH92. 



3. Tuning the Hadron Cocktail Packages 

To tune and validate a given hadron-cocktail package 
we can compare its output to three measured hadron 
distributions: 

• The pt distribution of tracks that stop in MuID 
Gap2 (counting from 0), with p z larger than a given 
minimum value. 

• The pt distribution of tracks that stop in MuID 
Gap3 (counting from 0), withp z larger than a given 
minimum value. 

• The vertex z distribution of reconstructed tracks, 
normalized to the collision- vertex distribution. 

Particles that stop in MuID Gap2 or Gap3 are those 
tracks for which no hit is found in the downstream gaps 
(Gap3 and/or Gap4). Figure [7] shows the longitudinal- 
momentum (p z ) distribution of tracks stopping in MuID 
Gap3 obtained using a given hadronic cocktail. Decay 
muons are characterized by a sharp peak, corresponding 
to electromagnetic energy loss in the absorber material. 
Note that the same peak would be obtained for muons 
from heavy-flavor decay. In contrast, hadrons are charac- 
terized by a broad shoulder that extends to much larger 
values of p z . For p z > p™ m (with p™ m « 3 GeV/c in 
this example) one obtains a clean hadron sample. The 
hadron-input pt distributions can then be tuned so that 
a good match between the number of stopped hadrons 
in the simulation and in real data is achieved in each pt 
bin. 

Figure [8] shows, for two muon-pT ranges, comparisons 
for real data and hadron-cocktail simulations of the z- 
vertex distributions of dN^/dz^BC tracks, which (a) are 
reconstructed in the north muon arm (located at positive 
z), (b) reach the MuID Gap4, and (c) are normalized by 
the event vertex distribution dN cvt / dz-QBC- The approx- 
imately linear dependence on zbbc is entirely due to the 
contribution of muons from light hadrons decaying before 
the muon-tracker front absorber. Muons from short-lived 
heavy-flavor hadrons have no measurable dependence on 
-Zbbc aud their contribution to the real-data sample is the 
source of the vertical offset between the hadron cocktail 
and the real-data distributions. Therefore, the hadron- 
cocktail package can be tuned by matching the slopes of 
these two distributions in each pt bin. The quality of 
this match is quantified by: 
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FIG. 7: (color online) Simulated p z distributions for particles 
that stop in MuID Gap3: (Black squares) all particles; (red 
triangles) stopped hadrons; (blue circles) decay muons. 
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FIG. 8: (color online) Vertex z distribution of tracks recon- 
structed in North (z > 0) MuID Gap4, for two transverse 
momentum bins. The real data (black closed circles) are 
compared to a given hadron-cocktail package (blue open dia- 
monds). The offset between data and the hadron cocktail is 
the contribution from heavy-flavor decays. 



c 0.4 
n 

a 

E 0.3 

1 0.2 
8 0.1 

C 

2 "° 

S-o.1 

n 

| -0.2 
|-0.3 
-0.4, 



South arm 



O FL105 
□ FL106 
A FL107 
<> GH91 
ft GH92 



V 



P T (GeV/c) 




P T (GeV/c) 



FIG. 9: (color online) Relative dispersion between the Nc 
yields obtained with the five hadron cocktails for the p + p 
analysis. Each hadron cocktail package is compared to the 
mean of the five packages for the north (top panel) and south 
(bottom panel) muon arm. 



Tuning of each hadron-cocktail package is achieved by 
iteratively selecting a set of pr-dependent weights (ap- 
plied to each track's thrown px) that simultaneously 
optimizes the agreement between data and simulation 
for the three distributions described above. Applying 
these weights to those simulated hadron tracks that reach 
MuID Gap4 determines the corresponding hadron contri- 
bution to the inclusive muon yield (Nc , Eq. [3]) . 

Figure |9] shows the relative dispersion between Nc val- 
ues obtained for the five different hadron cocktail pack- 
ages used for the p + p analysis, as a function of pr- For 
both muon arms, the largest differences exist between the 
gheisha and fluka cocktail packages for p T < 2 GeV/c, 
with a spread of about 20%. For p T > 3 GeV/c, most of 
the dispersion between the packages is due to increased 
statistical uncertainty in the data yields used to tune the 
hadron cocktail. 



Xg^(pt) = 



N ^(_AN 

i=l 



AN) 2 



07 



(8) 



where Nt,- ms is the number of zbbc bins; AN — 
dNi I dzQQc — dNc/dzBBC is the difference between the 
data and simulation for the i th zbbc bin; A7V is the av- 
erage difference over the entire zbbc range; Ci and cr mean 
are the statistical uncertainties of AN and AN, respec- 
tively. 



4- Systematic Uncertainties Associated with Individual 
Hadron Cocktail Packages 

There are two systematic uncertainties associated with 
the implementation of a given hadron cocktail package: 

csystPack the uncertainty associated with the imple- 
mentation of the hadron cocktail packages. It is 
comprised of two components: the uncertainty on 
the hadron cocktail input distributions and the so 
called MuID Gap3 to Gap4 matching uncertainty. 



11 



The uncertainty on the hadron cocktail input dis- 
tributions amounts up to 20% and is correlated be- 
tween the two arms. The uncertainty on the MuID 
Gap3 to Gap4 matching corresponds to tracks, in 
either real data or simulations, that get assigned 
an incorrect penetration depth, due to accidental 
addition of extra hits in the next MuID gap, or on 
the contrary, to detection inefficiencies. This un- 
certainty is evaluated using simulations. It is arm 
independent and amounts to 10%. These two con- 
tributions are uncorrelated and added in quadra- 
ture. 

cpackMismatch the uncertainty that characterizes, as a 
function of px, the ability of a given hadron cock- 
tail package to reproduce the measured distribu- 
tions described in the previous section. To eval- 
uate this uncertainty the cocktail is tuned three 
times, each time matching one of the three mea- 
sured hadron distributions perfectly. The disper- 
sion between the resulting background yields Nc 
obtained with these three different tunings, along 
with the central value for Nc obtained using the 
simultaneous tuning described above, is assigned 
to crpackMismatch- A different value is calculated for 
each muon arm, each px (and centrality) bin, and 
each of the five hadron cocktail packages. Math- 
ematical details of the calculation are outlined in 
Section IIII Gl Since the optimization is arm inde- 
pendent, this uncertainty is uncorrelated between 
the two muon arms. The magnitude of this un- 
certainty varies from 10 to 20% depending on the 
muon arm and the px bin. 

D. Other Background Sources 

In addition to the hadronic background, other back- 
ground sources include: 

• muons from heavy-flavor-resonance leptonic decay 
(e.g. Xc J/ip, ip' and the T family); 

• muons from Drell-Yan; 

• muons from light vector meson decay (ft, <fi an d <-o). 

These three sources contribute significantly less to the 
inclusive yields than the backgrounds from light hadrons. 
Monte Carlo simulations performed in the same manner 
as in [5| show that their contribution to the final heavy- 
flavor muon px spectrum is less than 5% in the px range 
used for this analysis and they have negligible impact 
with respect to the other sources of systematic uncer- 
tainties. 

E. Acceptance and Efficiency Corrections 

Acceptance and efficiency corrections, Ae, enter in the 
denominator of invariant yield measurements (Eq. [3]). 



TABLE II: Uncertainties in the acceptance and efficiency cor- 
rections. Individual components are added in quadrature to 



obtain the total value of ctae- 

Component Value 

OMuTr MuTr station data/MC 8% 

OMuiD MuID Gap4 efficiency uncertainty 4.5% 

o" run to run Run to run variation 2% 

o"p-scaie momentum scale 1.5% 

aAe Total 9.5% 



They are evaluated using simulated prompt single muons, 
propagated through the detector using the PHENIX 
GEANT3 simulation and reconstructed with the same 
analysis code and the same track quality cuts as for the 
real data analysis. These corrections account for the 
detector's geometrical acceptance and inefficiencies (for 
example, due to tripped high voltage channels or dead 
front-end electronic channels) . They also account for the 
muon triggers, reconstruction code and analysis cut inef- 
ficiencies. 

A reference run, representative of a given data taking 
period, is chosen to define the detector's response to par- 
ticles passing through it. This includes notably the list of 
inactive high-voltage and electronic channels. Remaining 
run-to-run variations with respect to this reference run 
are small due to the overall stability of the detector's 
performance, and are included in the systematic uncer- 
tainties (cr runtorun = 2%). 

A comparison between the hit distributions in the 
MuTr and the MuID obtained for the reference run in 
real data and simulations is used to assign an additional 
systematic error on our ability to reproduce the real de- 
tector's response in the simulations. Areas with unac- 
ceptable discrepancies are removed from both the simu- 
lations and the real data using fiducial cuts. Remaining 
discrepancies are accounted for with a 8% systematic un- 
certainty for the MuTr and 4.5% for the MuID. 

The hit multiplicity in the MuTr for Cu + Cu colli- 
sions is much higher than for p + p collisions and for the 
single muon simulations. To account for deterioration 
of the reconstruction efficiency in presence of such high 
multiplicity events, simulated single muon events are em- 
bedded into real data Cu + Cu events before running the 
reconstruction and evaluating the Ae correction. 

Another systematic uncertainty, cr p _ SC aic, is assigned 
to a possible systematic bias between the particle's re- 
constructed momentum and its real momentum. This 
uncertainty is estimated by comparing the measured J/ip 
invariant mass (using the dimuon invariant mass distri- 
bution) and its Particles Data Group (PDG) value. This 
uncertainty amounts to ~ 1.5%. 

Table [Til summarizes the acceptance and efficiency re- 
lated uncertainties, which sum quadratically to 9.5%. 
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TABLE III: Uncertainties in the single muon analysis. The 
individual components contribute to the final uncertainty as 
discussed in Section Till Gl 

Component Value 

cpackMismatch Package mismatch varies, ~10% 

osystPack Single package uncertainty varies, 10 - 20% 
avu Acceptance and efficiency 9.5% 
(Tbbc £ bbc 9.6% 

<jjv coll -Neon varies, 10 - 13% 



F. Systematic Uncertainties 

This section summarizes systematic uncertainties as- 
sociated with this analysis, most of which have been de- 
scribed in previous sections: 

• Systematic uncertainties associated with indi- 
vidual hadron cocktail packages, <7s ys tPack and 
cpackMismatch (Section |IIIC_4|) ; 

• Systematic uncertainty resulting from the dis- 
persion of the results obtained with the differ- 
ent hadron cocktail packages (introduced in Sec- 
tion |nrc21 mathematical details in Section [ill G[) ; 

• Systematic uncertainty on the acceptance and ef- 
ficiency correction factors, a At (Section IIIIEI and 
Table IJ); 

These systematic uncertainties are calculated indepen- 
dently for each arm, pt bin and centrality bin. 

The first three uncertainties listed above (first two 
items) are related to the hadronic background estimate 
and are combined to form a er mo dei systematic uncer- 
tainty, following a method described in Section [til Gl 

For invariant cross section measurements (in p + p col- 
lisions) and measurements of Raa one must add to the 
uncertainties above the systematic uncertainty on the 
p + p inelastic cross section seen by the minimum bias 
trigger, = 9.6%. For Raa measurements, one must 

also add the systematic uncertainty on the mean number 
of binary collisions (iV co ii) in each centrality bin, as pro- 
vided by the Glauber calculation used to determine this 
quantity. 

Table IIIII summarizes the systematic uncertainties in 
this analysis. 

G. Determination of the Central Value for 
Heavy-Flavor-Muon Production Yields and Raa 

This section details the procedure used to combine the 
results from multiple hadron cocktail packages to ob- 
tain the central values for the px spectra and -Raa and 
to propagate associated systematic uncertainties. This 
discussion includes the definition of <7p ac kMismatch and 



CModei- Throughout this section the variable Q is used 
to represent either the invariant yield or i?AA, for a given 
Pt and centrality bin; the procedure is the same for both, 
except where noted explicitly. 

1. For each px bin i, hadron cocktail package j, and 
package tuning k, we calculate the value Qij,k 
where: 

k—1 is the optimal tuning that best matches all 
three hadron distributions simultaneously (see 
Section |IIIC3j); 

k—1 is the tuning that best reproduces the pt dis- 
tribution of particles stopping in MuID Gap2; 

k—3 is the tuning that best reproduces the pt dis- 
tribution of particles stopping in MuID Gap3; 
and 

fc=4 represents the tuning that best reproduces 
the vertex z distribution of particles reaching 
MuID Gap4. 

The tuning k = 1 is used for the central value 
whereas the other tunings are used to establish the 
systematic uncertainty for a single hadron cocktail 
package due to its inability to completely describe 
measured hadron distributions. 

2. The package mismatch contribution to the uncer- 
tainty on the measurement Qij t k is estimated by 
the standard deviation between the four tunings, 
k: 



""PackMismatch.lJ — 4 / AQi<3,k (QiJ,k)) 2 (9) 
fc=l 

3. For each px bin i and package j, the associated 
total uncertainty cr^j is computed: 



a i,j ~ ^StatData.i + "StatPack,* ,j 

+ "SystPack.i + <7 PackMismatch,i ,j (1^) 
+ 

where the first two contributions, Og tatData , ; and 

"statPack i are * ne statistical uncertainties on the 
data and on the simulation and all other terms have 
already been introduced in previous sections. 

4. Using <7ij from step 3 we calculate the weighted 
mean of the Qij values obtained for the optimal 
tuning (k = 1) of the different packages, j, in each 
Pt bin, i: 

5 

(Qi) = /) w ijQi,j,k=l (11) 
J'=l 
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where 



IK, 
u ^ = - — — • 

3=1 



(12) 



5. The total uncertainty on the final measurement is 
the variance of the weighted mean: 



Var«Q i )) = XXX 



3 = 1 



2 X/ W *,3 W h™- a ' 



(13) 



common,? 



j<m 



where <7 C ommon,j is the part of the total uncertainty 
that is correlated between different packages: 



common. / 



CT StatData,i 
^SystPack,! + ° ' Ac, 



(14) 



6. For convenience, the total uncertainty Var^Q.;)) 1 / 2 
is split into a statistical uncertainty, a model- 
related systematic uncertainty and an acceptance 
and efficiency correction related systematic uncer- 
tainty: 



Var((Qi)) = 

^StatCombincd,? ~^ ^model,i + <?Ac,i ( 15 ) 

with: 



1 5 

^StatCombincd,! = ^StatData,; + g ^StatPack.j.j (16) 

3=1 



and (by construction): 

°modcU = Var((Qj)) — 0StatCombincd,i — a \c,i (1?) 

so that the final measurement, in a given muon 
arm, is written: 

(Qi) ± CTStatCombincd,i ± O"model,i ± CAe.i (18) 

7. The independent North and South measurements 
are combined using: 



(Qi) ^ ^ Wi,jQi,j i 
3=1 



(19) 



where i is the index of the pt bin, j the arm index 
and wij a weight calculated in the same manner as 
in Eq. 1121 using the following total uncertainty on 
the measurement Qi,f 

= "StatCombincd^j + ^model.i.j + a Ac,i,ji (20) 



which is identical to the expression of Eq. [T5J but 
explicitly includes the arm index, j. 

The total uncertainty on the arm-averaged Qi value 
is calculated in a manner similar to Eq. 1131 



Var«Q i » = X ■'■■,,<, 



3=1 



2 X W h3 W i,m cr , 



(21) 



arm common,? 



j<7n 



where er 2 rm common j is the systematic uncertainty 
common to both muon arms due to uncertainty on 
cocktail input. 

For convenience, this uncertainty is again split into 
a statistical contribution <7statCombined,i and a sys- 
tematic contribution <7s ys tCombined,i defined by: 

^SystCombined.i = Var((Qj)) — 0"statCombined,i ( 22 ) 

so that the final, arm-averaged, measurement of Qi 
is written: 



(Qi) ± (TSystCombined.i ± CStatCombi] 



(23) 



As already noted in Section IIIIF1 for invariant cross 
section measurements (in p + p collisions) and Raa mea- 
surements one must add the systematic uncertainty on 
the p + p inelastic cross section seen by the minimum bias 
trigger, cr^ BC in quadrature to the uncertainties above. 
For Raa measurements one must also add the system- 
atic uncertainty on the mean number of binary collisions 
-/V co n in each centrality bin. 



IV. RESULTS 

A. Heavy-Flavor Muon pr Distributions in p + p 
and Cu + Cu Collisions 

Figure [TU] shows the production cross section of neg- 
atively charged muons from decays of open-heavy-flavor 
mesons as a function of pt- Vertical bars correspond to 
statistical uncertainties. Boxes correspond to the sys- 
tematic uncertainties calculated following the steps de- 
scribed in Section MI Gl The measurements from both 
muon arms have been combined to reduce the overall un- 
certainty. Measured values for each pt bin are listed in 
the Appendix (Table EJ). 

Figure [TT] shows the invariant yield of negative 
muons from heavy-flavor mesons decay for all analyzed 
Cu + Cu centrality classes, compared to the invariant 
yield measured in p + p collisions. The solid lines cor- 
respond to a fit to the p + p data using the function 
A(l + (p T / B) 2 )~ 4 - 2 , similar to the one used in [35| . 
scaled by the average number of binary collisions iV" co n 
for each Cu -I- Cu centrality bin. For peripheral (40 — 
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FIG. 10: (color online) Production cross section of negative 
muons from heavy-flavor mesons decay as a function of pr in 
p + p collisions at \/s = 200 GeV. 
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FIG. 11: (color online) Invariant production yields of neg- 
ative muons from heavy-flavor-mesons decay as a function 
Pt in p + p collisions at i/s = 200 GeV (open squares) and 
in Cu + Cu collisions for three different centrality intervals 

(40 94%, 20 40% and 20%), scaled by powers 

of ten for clarity (filled circles). The solid line associated to 
each set of points corresponds to a fit to the p + p invariant 
yield distribution described in the text, scaled by the appro- 
priate number of binary collisions JV C oii when comparing to 
the Cu + Cu measurements. 



—94% centrality) and midcentral (20 40% centrality) 

Cu + Cu collisions, a reasonable agreement is observed 
between the measurement and the scaled fit to the p + p 
data, whereas for central collisions (0 20%), a sys- 
tematic difference is visible for high px muons (px > 
2GeV/c), and the measurements are below the scaled 



p + p fit, indicating a suppression of the heavy-flavor 
yields with respect to binary collision scaling, which is 
best quantified by measuring Raa (see Section HVC|) . 



B. Charm Cross Section, da C c/dy\ i y \ =1 65 in p + p 
Collisions 

The p + p heavy-flavor muon px distribution is used 
to estimate the charm differential production cross sec- 
tion, da C c/dy at forward rapidity ((y) = 1.65), as de- 
scribed in detail in reference [27] . The muon px spec- 
trum measured in p + p collisions spans from px = 1 
to 7 GeV/c. Estimation of the full charm charm cross 
section requires a theoretical calculation in order to ex- 
trapolate the measurement down to px — GeV/ c. A set 
of fixed-order-plus-next-to-leading-log (FONLL) [H, H3] 
calculations have been used in this analysis. 

The charm production cross section da cE /dy is derived 
from the heavy-flavor muon cross section using: 



da cS /dy 



1 



1 



da,, 



BR(c -> fi) C, 



l/D 



dy 



(24) 



where BR(c -+ D) is the total muon branching ratio of 
charm and is fixed to 0.103 in FONLL; Ci/o is a kine- 
matic correction factor, also provided by the FONLL cal- 
culation, which accounts for the difference in rapidity 
distributions between leptons and D mesons; dcr„- j dy 
is the total cross section for negative muons from heavy- 
flavor mesons decay, integrated over px and estimated by 
extrapolating our measurement down to px — GeV/c 
using the FONLL calculation. 



1. Extrapolation of the Data for px < 1.0 GeV/c 

Low px muons dominate the integrated heavy-flavor 
muon cross section due to the power-law like behavior of 
the px distribution (Fig. [TU|): according to the central 
value of the FONLL calculation, the integrated charm 
cross section for px^ > 1 GeV/c represents about 6% 
of the total. Additionally, the contribution of bottom 
quark decay to the heavy-flavor muon px distribution 
becomes increasingly relevant for px > 4 GeV/c, but has 
a negligible contribution to the integral and is ignored 
hereafter. 

The measured spectral shape matches the calculated 
shape. Therefore, extrapolation of the measured heavy- 
flavor muon px spectra down to px = GeV/c using 
FONLL is given by: 



dcr cE /dy| PHENIX = dcF c5 /dy\ 



FONLL ^FONLL ^ 



where q, fonll is a constant determined by fitting the 
central values of the FONLL px distribution to the data 
for px > 1 GeV/c. ft amounts to 3.75, and is used in 
determining the central point for PHENIX muons shown 
in Fig. HU 
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2. Systematic Uncertainties on dcr C c/dy\^ y y =1 65 

The total systematic uncertainty assigned to 
da C c/dy\/y\—i 65 is a combination of experimental 
and theoretical uncertainties, added in quadrature. The 
experimental systematic uncertainty on the integral 
above pr > l,GeV/c is determined by the appropriate 
quadrature sum of the systematic uncertainties on the 
individual px points. This uncertainty is up/down 
symmetric and is equal to 32 %. 

The theoretical uncertainty for da C c/dy\i y \_ x 65 origi- 
nates from the FONLL uncertainties. The variation in 
the FONLL calculation are determined by variation of 
the factorization scale, fiF, the renormalization scale, hr, 
and the charm quark mass. Other contributions, such 
as fragmentation and parton distribution functions are 
smaller and neglected in this analysis. 

The FONLL upper and lower bounds obtained by vary- 
ing the scales and the charm quark mass are treated as 
approximations for a one standard deviation systematic 
uncertainty. The ratio of the measured pt distributions 
for pt > 1 GeV/c to the upper and lower FONLL bounds 
are fit independently to determine the corresponding two 
normalization factors. The difference between these two 
normalization factors is then used as a theoretical un- 
certainty. This uncertainty is asymmetric and amounts 
to %. These FONLL systematic uncertainties are 
consistent with those of a previous study [27j , which ex- 
amined the different pr distributions obtained by vary- 
ing the FONLL parameters, 1.3 < M c [GeV/c}< 1.7, 
0.5 < ^iii/nriT < 2, 0.5 < /if/^t < 2, with tot rep- 
resenting transverse mass. The different predicted theo- 
retical px ranged within an envelope of ±35% relative to 
the central spectrum. 



3. Integrated Charm Production Cross Section at 
(y) = 1.65 in p + p collisions 

The integrated charm production cross section at for- 
ward rapidity ((y) = 1.65) obtained with this method 
is: 



p+p\/s = 200 GeV 
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FIG. 12: (color online) cc production cross section as a func- 
tion of rapidity in p+p collisions, measured using semileptonic 
decay to electrons (closed square) and to muons (closed cir- 
cle). 



C. Heavy-Flavor-Muon Raa in Cu + Cu Collisions 
as a Function of pr 

Figure [T3l shows Raa(pt) for muons from heavy- flavor 
decay in Cu + Cu collisions as a function of muon px 

for three centrality classes (40 94%, 20 40% and 

20%). As was the case for invariant yields and 

cross sections, the two independent measurements ob- 
tained with each muon arm are statistically combined, 
following the method discussed in Section MI Gl Verti- 
cal bars correspond to the statistical uncertainties; boxes 
centered on the data points correspond to point-to-point 
correlated uncertainties and the vertical gray band cen- 
tered on unity corresponds to the uncertainty on N co u, 
as listed in Table HI Also shown in the bottom panel of 
Fig. EH is a theoretical calculation from [3{||4(|, discussed 
in Section |Vj The measured values for each pt bin and 
each centrality class are listed in the Appendix (TablefVT). 



V. DISCUSSION AND CONCLUSIONS 



d<?cc/dy\ {y)=1 65 = 0.139 ±0.029 (stat) +°'°^ (syst) 

(26) 

This measurement is shown in Fig. 1121 together with the 
measurement performed by PHENIX at midrapidity [38[ , 
as well as the FONLL calculation and its uncertainty 
band, calculated as discussed in the previous section. 
The full circle, located at y — —1.65, corresponds to the 
combined measurement performed in both muon arms. 
The open circle, located at y — 1.65, corresponds to its 
mirror image. 



The measurement of open-heavy-flavor muon produc- 
tion in p + p collisions at ^fs — 200 GeV reported in 
this paper is a significant improvement over the previ- 
ous PHENIX published result [25T |. The transverse mo- 
mentum range of the present measurement is extended 
to pr = 7 GeV/c (compared to = 3 GeV/c in 
the previous analysis). The differential production cross 
section is integrated over pt to calculate a production 
cross section at forward rapidity of da cS / dy\^ =16b = 

0.139 ± 0.029 (stat) ^p^g (syst) mb. This cross section 
is compatible with a FONLL calculation within exper- 
imental and theoretical uncertainties. It is also com- 
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FIG. 13: (color online) Transverse momentum distribution 
of 7?aa for negative muons from heavy-flavor mesons decay 
in Cu + Cu collisions in the following centrality classes: 40- 
94% (top panel), 20-40% (middle panel) and 0-20% (bottom 
panel). Also shown in the bottom panel is a theoretical cal- 
culation from [39|, |40J , discussed in Section [V] 



FIG. 14: (color online) Comparison of Raa as a function of 
A'part between heavy-flavor muons reconstructed at forward 
rapidity (1.4 < y < 1.9) and p T > 2GeV in Cu + Cu colli- 
sions (red squares) and nonphotonic electrons reconstructed 
at midrapidity and pr > 3GeV in Au + Au collisions (blue 
circles). 



patible with expectations based on the corresponding 
midrapidity charm production cross section measured by 
PHENIX. 

Muons from heavy-flavor decay have also been mea- 
sured in Cu + Cu collisions at -Jsnn = 200GeV/c, in 
the" same rapidity and momentum range. This allows de- 
termination of the heavy-flavor-muon Raa as a function 

of p T in three centrality classes, 40 94%, 20 40% 

and 20%. As shown in Fig. Q21 no suppression is 

observed across most of the transverse momentum range 

for muon yields measured in peripheral (40 94%) and 

midcentral (20 40%) Cu + Cu collisions compared to 

-/V co n-scalcd p + p collisions. On the contrary, open heavy 
flavor production is significantly suppressed for central 

Cu + Cu collisions (0 20%), with the largest effect 

observed for pr > 2 GeV/c. Interestingly, as demon- 
strated in Fig. [T31 the level of suppression for these higher 
Pt heavy- flavor muons (the last red point on right) is 
comparable to the level of suppression observed for high 
Pt nonphotonic electrons measured at midrapidity in the 
most central Au + Au collisions (the last blue point on 
right). One expects the Bjorken energy density of the 
matter produced in the midrapidity region in the most 
central Au + Au collisions to be at least twice as large 
as that of the matter produced in the forward rapidity 
region in most central Cu + Cu collisions [26J, |4l| . There- 
fore the large suppression observed in Cu + Cu collisions 
suggests significant (cold) nuclear effects at forward ra- 
pidity in addition to effects due to strongly interacting 
partonic matter. 

As shown in the bottom panel of Fig. [T3l the suppres- 
sion of open-heavy-flavor muon production for central 
Cu + Cu collisions is consistent with a recent theoretical 
calculation performed at the same rapidity (y = 1.65) for 
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p T > 2.5 GeV/c [H,|4Ci|. This calculation includes effects 
of heavy-quark energy loss (both elastic and inelastic) 
and in-medium heavy meson dissociation. Additionally, 
the calculation accounts for cold nuclear matter effects 
relevant for open heavy flavor production [421 ] . namely 
shadowing (nuclear modification of the parton distribu- 
tion functions of the nucleon) and initial state energy 
loss due to multiple scattering of incoming partons be- 
fore they interact to form the cc pair. 

New PHENIX inner silicon vertex detectors will 
greatly improve heavy flavor production measurements 
and allow separation of charm and bottom contribu- 
tions mm. 
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APPENDIX: DATA TABLES 

Table IIVI gives the differential invariant cross section 
of muons from heavy-flavor decay in y/s = 200 GeV 
p + p collisions and corresponds to Fig.[TUJ Table IVl gives 



i?AA of muons from heavy-flavor mesons decay for the dif- 
ferent centrality classes of snn = 200 GeV Cu + Cu col- 
lisions and corresponds to Fig. (T3J 

TABLE IV: Differential-invariant cross section of negative 
muons from heavy-flavor mesons decay for 200 GeV p + p 
collisions at midrapidity. 

Pt (GeV/c) l/27rpT d 2 a/dpTdr) (mb) stat error syst error 



1.12 


3.64o-04 


1.55e-05 


1.23c-04 


1.36 


1.19c-04 


2.85o-06 


3.39e-05 


1.61 


4.57c-05 


1.15o-06 


1.25e-05 


1.86 


1.92e-05 


5.39e-07 


5.18c-06 


2.11 


8.31c-06 


3.02e-07 


2.27c-06 


2.36 


3.52e-06 


1.54e-07 


1.18c-06 


2.61 


1.67c-06 


9.21c-08 


5.74c-07 


2.86 


9.12c-07 


6.04e-08 


3.18c-07 


3.21 


3.83c-07 


2.28e-08 


1.22e-07 


3.72 


1.41c-07 


1.24e-08 


4.15c-08 


4.38 


3.34c-08 


3.49c-09 


1.12c-08 


5.65 


2.99e-09 


1.09o-09 


1.31e-09 



TABLE V: Nuclear-modification factor, Raa, of negative 
muons from heavy-flavor mesons decay as a function of pr 
for the specified centrality classes of Cu + Cu collisions at 
= 200 GeV. 



Centrality 


p T (GeV/c) 


Raa 


stat error 


syst error 




1.13 


6.93o-01 


3.98c-02 


1.87e-01 




1.38 


5.41e-01 


3.49e-02 


1.87o-01 




1.63 


6.57o-01 


5.32e-02 


2.20e-01 


0-20% 


1.875 


6.26e-01 


6.74e-02 


2.28e-01 




2.25 


4.54c-01 


6.90e-02 


1.50e-01 




2.75 


3.61c-01 


1.09e-01 


1.46e-01 




3.5 


3.95c-01 


1.46c-01 


2.00e-01 




1.13 


1.03c+00 


5.59e-02 


2.66e-01 




1.38 


9.32c-01 


4.63e-02 


2.46c-01 




1.63 


l.lle+00 


6.95c-02 


3.72c-01 


20-40% 


1.875 


1.34e+00 


1.08e-01 


4.59e-01 




2.25 


1.15c+00 


1.06e-01 


3.18c-01 




2.75 


8.14e-01 


1.40c-01 


2.80e-01 




3.5 


4.42o-01 


2.03e-01 


2.96e-01 




1.13 


1.36c+00 


7.27c-02 


3.38c-01 




1.38 


1.28c+00 


6.26e-02 


3.21e-01 




1.63 


1.16c+00 


8.87c-02 


3.08o-01 


40-94% 


1.875 


1.16c+00 


1.37e-01 


3.30c-01 




2.25 


8.64c-01 


1.43e-01 


2.73c-01 




2.75 


6.94e-01 


1.92e-01 


3.30c-01 




3.5 


8.09c-01 


2.80e-01 


3.47e-01 
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